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Abstract 



A two-gluon-exchange model incorporating perturbative and non-perturbative effects 
is presented for p electroproduction which provides an excellent description of all cur- 
Q_ij rent data. This is then used to calculate the contribution from the p to deeply virtual 

Compton scattering via the vector-meson-dominance transition p — > 7. This is found to 
be sufficiently large to provide a significant contribution through interference with the 
perturbative QCD term. 
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1 Introduction 

There is considerable evidence in 7*p reactions that the nominally perturbative regime can 
be strongly influenced by non-perturbative effects. This is an obvious feature of recent dipole 
models of deep inelastic scattering |], ]2j, ]3|], where for transverse photons especially the contri- 
bution from large (non-perturbative) dipoles extends to significantly large values of Q 2 . The 
penetration of non-perturbative physics into the perturbative regime is even more explicit in 
generalised vector dominance models ]!] or in two-component models ||, || |7|, [8], || |1(| which 
combine "soft" (non-perturbative) and "hard" (perturbative) contributions. Typically the soft 
contribution comprises the normal reggeon and soft pomeron exchanges, the latter with an 
intercept of ~ 1.08. The hard contribution may be a second pomeron, the hard pomeron, with 
an intercept of ~ 1.44 || 0, ||, or be based explicitly on perturbative QCD || |], |10|] . 

A good illustration of the two-component approach is provided by exclusive p electroproduction, 



7*p — > pp. The high-energy data |TT], [12] indicate that this approach is appropriate, as the 
effective pomeron intercept increases from the canonical hadronic value of ~ 1.08 for real 
photons to perhaps as large as ~ 1.19 at Q 2 = 20GeV 2 . In section 2 we present a two- 
component model for p electroproduction which successfully describes all current data. We 
base our calculations on two-gluon-exchange models of the pomeron. For the non-perturbative 



contribution we use the model of Diehl [13| and for the perturbative contribution that of Martin, 
Ryskin and Teubner [Q. The procedure follows the suggestion of |14| by calculating the light 
quark anti-quark pair uu and dd production process 7*p — > qqp with the invariant mass of 
the qq-p&ir Mx afterwards integrated over the mass interval of the p. This approach has the 
benefit of avoiding vector meson wave function complications, which can be serious ||15| , and 
allows one to concentrate on the production dynamics. 

A related topic is that of Deeply Virtual Compton Scattering (DVCS) on protons, 7*p — > jp, 
which is seen as an important reaction for the study of diffraction in QCD. In the standard 
QCD approach the amplitude is described by skewed parton distributions |16|] corresponding to 
operator products evaluated between protons of unequal momenta. These are generalizations of 
the familiar parton distribution of deep inelastic scattering, and like them satisfy perturbative 



evolution equations |T7[ which enable them to be evaluated at all Q in terms of an assumed 



input at some appropriate Q 2 = Q 2 ,. Preliminary data pi have been presented which are 



consistent with QCD predictions |TJ|, subject to two uncertainties 



The first is that the theoretical predictions refer to zero momentum transfer t = and to 
compare with experiment, one must integrate over t. This is done by assuming an exponential 
dependence exp(— bt) and estimating the unmeasured slope parameter b. The considerable 
uncertainty^ in b leads to a corresponding uncertainty in the normalization of the predictions. 

Secondly, it is necessary to specify the input skewed parton distributions at the reference 
Qo- In HH these are obtained by estimating their ratio to "ordinary" parton distributions 



tin analysing their preliminary data, HI |18 assume 7 < b < 10 GeV 



at Qq = 2.5 GeV 2 using arguments based on the aligned jet model pOfl, which in practice is 



almost identical to the simplest diagonal generalized vector meson dominance model for the 
soft pomeron term |2T[. The resulting ratio is of the order of 2, leading to a factor of order 
four in the predicted DVCS cross-sections. While this provides a reasonable first estimate, it 
is clearly subject to uncertainties which will become important when more accurate data are 
available. 

Here we note that there are direct "hadronic" contributions to DVCS via the vector-meson- 
dominance mechanism 7*p — > Vp, V — > 7. One particular vector- meson contribution to DVCS, 
namely that of the p, can be calculated with reasonable precision using the results of Section 2. 
This is done in section 3, where we show that the results provide useful constraints on models 
used to estimate the skewed parton distributions at the reference Q^. 



2 Rho electroproduction 



According to the factorisation theorem [2^] the exclusive vector meson production processes 
can be factored into three parts: the fluctuation of the (virtual) photon into a gg-pair; the 
interaction of the gg-pair with the proton; and the formation of the vector meson from the 
gg-pair. Similarly, the upper part of either diagram in figure [1], containing the process 7* — > gg 
can be considered separately from the rest of the diagram. Apart from the couplings a s , 
the upper parts of the non-perturbative and perturbative diagrams are identical. The lower 
parts, that is the parts containing the gluons and the proton, are described differently in each 
approach. The model of Diehl JT3] follows Landshoff and Nachtmann (LN) p5| |. The gluons 
are assumed not to interact with each other and a non-perturbative gluon propagator fl3| is 
used: 



D-np \ ) ~ A/np 
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4. The normalisation J\f np is determined from the condition 
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The phenomenological parameters /?o, which describes the effective coupling of the pomeron to 
the proton, and fio are determined from the total pp and pp cross section data and from deep 
inelastic scattering: j3 ps 2.0 GeV -1 and /i ~ 1.1 GeV [p4] . For the non-perturbative couplings 
of the gluons to the quarks forming the p a value a s « 1 is taken. 

It has been argued |23j that the diagrams in which the non-perturbative gluons couple to 
different valence quarks in the proton are suppressed and therefore can be disregarded. Only 
the diagrams where both gluons couple to the same valence quark are calculated. Each of the 
three valence quarks is incorporated into the proton according to the Dirac form factor of the 
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Figure 1: One of the four diagrams corresponding to the non-perturbative model [13|, |23[ (upper 
diagram), at t = A 2 , and the perturbative description |1 | (lower diagram), at t — 0. The other 
three diagrams differ in the way gluon lines are attached to the quarks in the qq-loop. The 
minus sign in the four-momenta indicates an antiparticle. The off-shell quark and the cut along 
which the quark lines are put on-shell in the LN model are indicated by the cross and the 
dashed line respectively. 



proton F lp (t), where t = A 2 and the four-momenta of the particles are as depicted in figure [J]. 
The energy dependence due to the soft pomeron comes via a factor x P aF in the amplitude, 
with x P = {M\ + Q 2 - t)/(W 2 + Q 2 - m 2 roton ) and a P (t) the soft pomeron trajectory pSfl . In 
principle we can calculate the ^-dependence of the soft-pomeron contribution from this, but as 
this cannot be done for the perturbative contribution, we calculate only at t — and use the 
experimental slope to give the integrated cross section. 



Following the argument of |13| , the coupling a s at three of the four vertices is taken at a non- 



„(o) 



perturbative scale, i.e. a s is used, while for the vertex where the gluon couples to the off-shell 
quark it is taken at a perturbative scale A 2 = 
scale for the whole upper part of the diagram. 



quark it is taken at a perturbative scale A 2 = (l 2 + m 2 ) (Q 2 + M^) /M|, which is a typical 



In the perturbative approach by Martin et al. [14]] the pomeron is modelled as a pair of per- 
turbative gluons with symmetric momenta. The perturbative gluon propagator V p (k 2 ) = 1/k 2 
is used. In principle the gluon flux can be obtained from the unintegrated gluon density 
f(x P , \k 2 \), which gives the probability of finding a t-channel gluon with the momentum squared 
\k 2 \ in the proton. However, a special treatment of the infrared region is required as the un- 
integrated gluon density f(x P , \k 2 \) is theoretically undefined as \k 2 \ — > and numerically 
unavailable below some value of \k 2 \ = Q 2 ,, which varies with the parton distribution chosen 
and usually is in the region from 0.2 to a few GeV 2 . The linear approximation as suggested in 



|i~4| is used to account for the contribution to the integral from the \k 2 \ < Qq region. A num- 
ber of gluon distributions were tried, including MRS(Rl), MRS(R2), GRV94HO, GRV94LO 



and others using the PDFLIB program libraries |26] for numerical calculations. However, the 
one which gave the best energy dependence within the model is CTEQ4LQ p7|, so results are 
presented only for that choice. 

The derivations of both models can be expressed in a common kinematical framework, taking 
into account the on-shell conditions along the cut line in figure [I], which result in the gluon 
momenta being predomimently transverse with respect to the 7*p axis, \k 2 \ m k t . Here and 
subsequently transverse two-vectors are shown in bold. Integrating over the azimuthal angles, 
one obtains a common structure for both models 103 at t = 0: 
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where Q = m 2 + Q 2 y£ 2 + m 2 ) /Ml and the symbol V denotes the model-dependent parts 
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Here f(x P , k 2 ) is related to the gluon distribution g(x F , Q 2 ) by 



r-Q 2 dk 2 
x P g[x P ,Q ) = / i_2 J\ x pi kt, 



k 



(8) 



with the inverse 



f[Xp, k t ) — fe ; 



o[Xp, g[Xp, k, 



dk 
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Such a common structure should be present since both models describe the same physical 
process. The common parts originate mainly from the kinematics of the process. The difference 
is contained in equations (6) and (7), arising from the different physical interpretations of the 
internal dynamics of the process in the two models. 

The relation between the models can formally be written as a replacement V np < — > V p . There 
are no tJqZ, couplings for the two bottom vertices in the perturbative case since the gluons are 



considered as part of the proton and described by fixp, k 2 j . A minor difference is the different 
argument of a s in both models and the fact that there is no need for the linear approximation 
in the nonperturbative model since the integration over the gluon momentum can be performed 
down to zero. The expressions above are given in a general form but can be further simplified 
for the light quarks assuming m q = 0. 

The forward differential cross section da/dt\t=o can be related to the total cross section a(W, Q 2 ) 
using the experimentally measured forward diffractive slope b p (Q 2 ) assuming an exponential 
^-dependence of da/dt: 



a 



7*p->pp 
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(10) 



t=o 



The slope parameter b p (Q 2 ) varied from 7 at the smallest Q 2 to 4 at the highest Q 2 . The 
M^-integration limits Mf = (0.6 GeV) and M 2 = (1.05 GeV) were chosen to span the p- 



region, following |14j . The polarisation of the photon beam e expt is a known characteristic of 
the experiment. For HERA e expt ~ 1 while for fixed-target experiments it varies significantly, 
depending on the energy and photon virtuality. It is generally in the range of 0.5 to 0.9 and 
this is taken into account when comparing our results with the data. 

Neither model by itself can describe all the observed features of the p electroproduction data 
simultaneously: that is, the absolute value of the total cross section a(W, Q 2 ); its dependencies 
on Q 2 and W; and the variation of the longitudinal to transverse cross section ratio R = a L jo T 
with Q 2 . The non-perturbative approach gives an energy dependence which is too flat at the 
higher values of Q 2 due to the soft pomeron intercept while the energy dependence of the 
perturbative approach, coming from the energy dependence of the gluon density, is clearly too 
steep at all Q 2 . The perturbative model does not replicate the Q 2 -dependence of the NMC data, 
and the non-perturbative model gives a longitudinal to transverse ratio which is somewhat low|. 

The results of summing the perturbative and non-perturbative production amplitudes 



V 



V 

I np 



v„ 



(11) 



with V np and V p as given by @ and (0), are shown in figures 0, [| and f| together with the data. 
It is clear that the two-component model gives excellent agreement. The effect of the value of 
Ql is shown in figure 0. The effect is consistent both at high and low energy and allows us to 
fix Ql w 0.9 GeV 2 , which is a reasonable value. 



*For further details, see p6 
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Figure 2: The dependence of total cross sec- Figure 3: The Q 2 -dependence of the longitu- 
tion on Ql compared with high energy data dinal to transverse ratio compared with high 



from HI |LJ and ZEUS 0, gj |2J, and low- energy data from HI fig |4) |§ and ZEUS 
energy data from NMC |£|. [|TT|, [32|, |36]], and low energy data from CHIO 
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Figure 4: VU-dependence compared with data from HI [1^], ZEUS JTT| and NMC [ 33fl . The use 
of the open and filled markers is only to indicate to which value of Q 2 the data points belong. 
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P 

Figure 5: The hadronic contribution to DVCS arising from a virtual p in the final state. 

3 Deeply virtual Compton scattering (DVCS) 

We now use the results of the previous section to estimate the contribution to DVCS from the 
mechanism of figure [|; and then comment on its implications for the estimation of the skewed 
parton distributions at the reference Q\ = 2.6 GeV 2 . 

3.1 Estimating the p contribution 

Assuming s-channel helicity conservation, only transverse photons contribute to the DVCS 



cross section. The relation between the DVCS and the ep — > cyp cross sections is [[29 

^ O pn — >prvr> ex Prri '' 



' ep^ejp 



dWdQ 2 n Q2 (W2 + Q 2 _ ^^ 



l + (l-y) 2 a^ w , (12) 



with y = [W 2 + Q 2 — m 2 roton ) / (s — rn 2 ro ton)^ V$ * s the centre-of-mass energy of the ep system 
and W is that of the 7*p system. The contribution from the mechanism of figure |5| is then 
given by 

Alio/. 
o P {l*P -* IP) « —j- <?{l*P -* PP) • (13) 

The coupling e/-y p of the p to the photon is directly related to the width of the decay p — > e + e~. 
In the narrow-width approximation 



4vr 3 T 



p^e T e" 



lp m p a 2 



0.494 ±0.023 (14) 



where the experimental values of m p , T p and the branching ratios have been used |30|| . 

The values obtained from (|T2"D, ([13|) and (|I4"D for the differential cross section da p /dQ 2 , aver- 
aged over each Q 2 -bin and integrated over 30 GeV < W < 120 GeV are shown in Table 1, 



together with the preliminary experimental values |T8[ for da/dQ with the Bethe-Heitler term 



subtracted off. As calculated, the amplitude in the model is purely imaginary^. We also give 
the estimate 



RAQ 



Ap (Yp -> IP) 



da p /dQ 2 
da/dQ 2 



1/2 



(15) 



of the amplitude of figure || to the total non-Bethe-Heitler amplitude. As can be seen, this 
reduces from of order 20 % in the lowest Q 2 bin to of order 10 % in the highest. Although 
small, this is not a negligible contribution to the total amplitude for DVCS, since it will interfere 
constructively with the remaining dominant contributions assuming they are mainly imaginary. 
It is also worth noting that even for the lowest Q 2 bin, A p (j*p — >■ 7p), which one might expect to 
be predominantly "soft", contains a significant contribution, about 40 %, from the perturbative 
term. 



30GeV <W< 120 GeV 



Q 2 bin 



GeV' 



2.0 to 4.0 

4.0 to 6.5 

6.5 to 11.0 

11.0 to 20.0 



da p /dQ 2 
pb/GeV 2 



1.576 

0.213 

0.0352 

0.00498 



da expt /dQ 2 
pb/GeV 2 



47 



+12 
10 



2.10 
0.35 



0.51 

0.64 

+0.17 

0.14 



Ra 



1 O+0.02 
io -0.02 





0.18 
0.13 
0.12 



0.04 

0.02 
+0.02 

0.02 
+0.02 

0.03 



Table 1: Comparison of the preliminary experimental values of da/dQ 2 after subtraction of 
the Bethe-Heitler contribution, integrated over 30 GeV < W < 120 GeV, with the hadronic 
contribution of figure [5]. The corresponding ratio of amplitudes Rj, is also given, assuming the 
amplitudes are pure imaginary. 



3.2 Estimating the input distributions 



In |I9| the input skewed parton distributions at the reference Qq = 2.6 GeV were estimated by 
relating them to "ordinary" parton distributions by using arguments based on the aligned jet 
model ]2U| . In practice this is almost identical to the application of generalised vector meson 
dominance (GVD) in its diagonal form [21 1. Both models assume that for the scattering of 
virtual photons the amplitude at t = is of the form 



m 2 p(m 2 



l dm (Q 2 + m 2 ) 2 



(16) 



111) 



Im„4( 7 *iV^7*A0 t= o = — I 
§There is no interference term, since the data are integrated over the azimuthal angle 4> r . (See eqn.(40) of 



where p(m 2 ) is taken to be energy independent, corresponding to soft Pomeron behaviour with 
intercept a P = 1. For the aligned jet model p(m 2 ) is given by 



P(m 2 ) = a^ M (m 2 )R e+e ~(m 2 ) 3 -^. (17) 



In (jTjJ the product a r ^f M (k^)R e e (m 2 ) is assumed to be a constant. For diagonal vector meson 
dominance p(m 2 ) is given by 

p{m 2 ) = a™(m 2 )R e+e -(m 2 ), (18) 

and the identification of the two approximations is completed by the usual diagonal GVD 
assumption that <jy°p(m 2 ) ~ 1/m 2 for a P = 1 and R e e (m 2 ) is constant. 

In the case of DVCS the imaginary part of the amplitude for t = is obtained from (|TJ]) by 
replacing one of the propagators with 1/m 2 . Since m 2 p(m 2 ) in (16) and in its DVCS equivalent 



is taken to be constant in the aligned jet/GVD model, the integrals are trivial, giving the result 
lmA{j*N->jN) t=0 1 



Im.4(7*iV ->7*iV)i = o "~ Q 



{m 2 + Q 2 ) ln(l + Q 2 /m 2 ) . (19) 



For a reasonable choice of the lower limit m 2 , in (|H3|), typically 0.4 to 0.6 GeV 2 , the ratio ~ 2 
for Q 2 « 2.5 GeV 2 . So knowing F 2 (x, Q 2 ) this gives the input to the DVCS evolution equations 



At this point we notice that the same aligned jet/GVD model also implies values for the 
p contribution to DVCS. Attributing the low mass contribution to the p-meson, one easily 
obtains 

ImA p ( 7 *iV->7A0 f =o ln(l + Q 2 /M 2 ) - ln(l + Q 2 /M 2 ) 

AW ' ImA( 7 *iV^7A0 <= o ' ln(l + Q 2 /M 2 ) { ] 

where we have chosen m 2 , = M 2 = (0.6 GeV) 2 and M^ = (1.05 GeV) 2 for consistency with 



our treatment of the p in Section 2 (cf. equation [TTJ). This gives values for Rj± of 0.43 at 
Q 2 = Ql = 2.6 GeV 2 reducing to 0.28 at Q 2 = 15 GeV 2 . These values are considerably larger 
than the much more reliable estimates of section 3.1. 

This discrepancy is not surprising given the extreme simplicity of the aligned jet/GVD model 
used. More elaborate versions can be constructed and will be required in the future to estimate 
more accurately the skewed parton distributions at the reference Qfc. Such models will need to 
take account of both "hard" and "soft" diffraction, and in constructing them the requirements 
that they be compatible with the p contributions of Section 3.1, as well as the structure function 
data for < Q 2 < Ql, will be essential constraints. 
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